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ABSTRACT 

A  new,  fully  conserved  minimal  adjustment  scheme  with  temperature  and  salinity  (T,  S)  coherency  is  presented 
for  eliminating  false  static  instability  generated  from  analyzing  and  assimilating  stable  ocean  (T,  S)  profiles  data, 
that  is,  from  generalized  averaging  over  purely  observed  data  (data  analysis)  or  over  modeled/observed  data 
(data  assimilation).  This  approach  consists  of  a  variational  method  with  (a)  fully  (heat,  salt,  and  potential  energy) 
conserved  conditions,  (b)  minimal  adjustment,  and  (c)  ( T,  S)  coherency.  Comparison  with  three  existing  schemes 
(minimal  adjustment,  conserved  minimal  adjustment,  and  convective  adjustment)  using  observational  profiles 
and  a  simple  one-dimensional  ocean  mixed  layer  model  shows  the  superiority  of  this  new  scheme. 


1.  Introduction 

Analysis  and  assimilation  of  ocean  temperature  and 
salinity  ( T ,  S )  profile  data  are  a  generalized  average  over 
purely  observed  data  (data  analysis)  or  over  modeled- 
observed  data  (data  assimilation)  (Chu  et  al.  2004).  For 
example,  data  assimilation  is  intended  to  blend  the  mod¬ 
eled  variable  (xm)  with  observational  data  (y0)  (Sun  1999): 

X«  =  Xm+W-[yo-tf(Xm)]’ 

where  x„  is  the  assimilated  variable,  H  is  an  operator  that 
provides  the  model’s  estimate  at  the  observational 
points,  and  W  is  the  weight  matrix  that  is  determined  by 
various  data  assimilation  schemes  such  as  optimal  in¬ 
terpolation  (e.g.,  Lozano  et  al.  1996),  Kalman  filtering 
(e.g.,  Galanis  et  al.  2006),  and  variational  methods  (e.g., 
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Tang  and  Kleeman  2004).  Due  to  the  nonuniform  vertical 
distribution  of  the  observational  profile  data  and  the  high 
nonlinearity  of  the  equation  of  state  of  seawater,  a  false 
static  instability  may  be  generated.  For  example,  10-day 
Jet  Propulsion  Laboratory  (JPL)  Estimating  the  Circula¬ 
tion  and  Climate  of  the  Ocean  (ECCO)  ( 7’,  S)  fields 
centered  on  31  December  2008  (download  on  19  February 
2009  from  the  JPL  web  site:  http://ecco.jpl.nasa.gov/ 
external/)  show  that  a  considerable  portion  (11.6%)  of 
the  profiles  are  statically  unstable  (Chu  and  Fan  2010a). 

For  discrete  samples  (Tk,  Sk)  at  depth  zk  (negative 
value),  k  =  1,  2,  . . .  ,  K  (k  increasing  downward),  the 
density  difference  between  two  adjacent  levels, 

Ek  =  P(Sk+ 1>  Tk+l’zk)  -  p($k’  Tk’zk^ 
k  =  1,2,  ...  ,K  —  1 ,  (1) 

is  taken  after  one  is  adiabatically  displaced  to  the  depth 
of  the  other.  Here,  Ek  is  defined  as  the  static  stability 
(Lynn  and  Reid  1968),  p(Sk+\,  Tk+i,Zk)  is  the  local  po¬ 
tential  density  at  depth  {zk+ 1)  that  the  water  parcel 
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would  acquire  if  adiabatically  brought  to  a  reference 
depth  ( Zk ),  and  p(Sk,  Tk,Zk)  is  the  in  situ  density  to  the 
depth  of  the  upper  of  the  two  adjacent  levels  ( Zk )•  The 
density  inversion  is  defined  by  the  occurrence  of  negative 
values  of  Ek.  The  minimum  static  stability  is  represented 
by  Ek  =  0.  It  is  not  always  possible  to  reach  zero  exactly 
due  to  the  precision  limitations  of  the  temperature  and 
salinity  values  used  (Locarnini  et  al.  2006).  As  a  result,  the 
minimum  value  for  the  static  stability  is  given  by 

Ek^E  mm’  *  =  1.2,...,*,  (2) 

where  Emin  is  the  reference  value  for  the  minimum  static 
stability,  which  is  a  user-defined,  small  positive  value.  If 
static  instability  occurs  in  observed  or  averaged  hydro- 
graphic  casts  [i.e.,  (2)  is  not  satisfied],  it  needs  to  be 
adjusted  (Chu  and  Fan  2010a). 

Let  an  unstable  pair  of  hydrographic  casts  [(7*.  Sk), 
k  =  1,  2,  . . .  ,  7T]  be  represented  by  a  2A-dimensional 
vector: 

x  =  (TvT2,  ...  ,Tk,SvS2,  ...  ,SK ).  (3a) 

An  adjustment  process  is  used  to  change  the  vector  x 
into  stable  casts: 

x  +  Ax  =  (T1  +  A Tv  T2  +  A T2,  ...  ,Tk  +  A T K, 

Aj  +  AAj ,  S2  +  AA2,  . . .  ,  A^  +  AA^) .  (3b) 

Ideally,  three  principles  should  be  considered  to  de¬ 
termine  Ax:  1)  conservation  of  heat,  salt,  and  potential 
energy;  2)  minimal  adjustment,  that  is,  minimizing  ||Ax|| 
subject  to  stabilizing  the  casts  x  +  Ax;  and  3)  preserva¬ 
tion  of  ( T ,  A)  coherency.  However,  these  principles  are 
not  fully  used  in  the  existing  three  types  of  schemes:  (a) 
convective  adjustment  (CA),  (b)  minimal  adjustment 
(MA),  and  (c)  conserved  minimal  adjustment  (CMA). 

The  CA  schemes,  often  used  in  ocean  modeling,  are 
based  on  the  same  original  idea  (e.g.,  Bryan  1969;  Yin  and 
Sarachik  1994):  whenever  a  water  column  is  statically 
unstable,  temperature  and  salinity  are  vertically  adjusted 
to  make  the  water  column  neutrally  stable,  with  heat  and 
salt  conserved  in  the  process.  The  adjustment  takes  an 
iterative  approach.  The  iteration  continues  between  all 
adjacent  levels  until  the  static  instability  is  removed  in 
the  entire  water  column.  Since  the  constraints  are  only 
the  conservation  of  heat  and  salt,  usually  the  adjustment 
|| Ax||  is  relatively  large  (Chu  and  Fan  2010b). 

The  MA  scheme  originally  proposed  by  Jackett  and 
McDougall  (1995)  is  employed  to  solve  the  problem: 

Minimize 1 1 Ax| | 2  =  AxTIAx  subject  to 

A  ■  (x  +  Ax)  >  E2mm ,  (4) 


where  I  is  the  identity  matrix.  The  finite-difference  ap¬ 
proximation  of  stability  Ek  becomes  the  inner  product  of 
the  matrix  A  and  the  profile  vector  x  +  Ax.  The  matrix  A 
depends  on  the  solution  Ax  to  the  minimization  problem 
(4),  implying  that  the  constraints  in  (4)  are  nonlinear. 

Generally,  the  norm  AxTIAx  in  (4)  is  not  desirable 
because  (a)  it  is  not  invariant  under  linear  transforma¬ 
tions  and  (b)  it  ignores  the  sources  of  the  errors  in  tem¬ 
perature  and  salinity.  A  quadratic  form,  AxTCAx,  with  C 
the  inverse  of  an  error  covariance  matrix  is  likely  to 
be  a  better  norm  since  it  will  yield  a  solution  in  which 
changes  of  trusted  values  tend  to  be  smaller  in  size  than 
the  changes  of  unreliable  values.  Furthermore,  the  qua¬ 
dratic  norm  is  a  natural  nondimensional  measure,  which 
can  easily  characterize  the  changes  in  temperature  and 
salinity.  However,  the  error  covariance  matrix  C  needs  to 
be  calculated  before  using  this  quadratic  norm. 

Usually,  an  iteration  method  is  used  in  the  MA  method. 
Based  on  the  MA  scheme  with  some  modifications, 
Locarnini  et  al.  (2006)  proposed  a  local  iterative  separated 
method  to  minimally  adjust  unstable  temperature  and  sa¬ 
linity  profiles.  Obviously,  the  MA  scheme  does  not  require 
conservation  of  heat,  salt,  and  potential  energy,  which  may 
lead  to  errors  in  estimating  the  oceanic  impact  on  global 
climate  change  patterns  (Chu  and  Fan  2010a)  since  oceanic 
heat  and  salt  transports  play  an  important  role. 

The  CMA  scheme  was  developed  by  Chu  and  Fan 
(2010a)  to  solve  (4)  by  requiring  conservation  of  heat  and 
salt.  A  set  of  well-posed  combined  (K+ 1)  linear  and  ( K —  1) 
nonlinear  algebraic  equations  for  { A  7T- AA;t,  k  =  1,2,  . . . , 
l<\  is  established  and  solved  using  Newton’s  method.  The 
CMA  scheme  has  the  three  features:  1)  conservation  of 
heat  and  salt,  2)  removal  of  static  instabilities  with  small 
( T ,  A)  adjustments,  and  3)  an  analytical  form.  However, 
it  does  not  conserve  the  potential  energy  and  preserve 
the  T-S  coherency  since  the  salinity  adjustment  is  as¬ 
sumed  to  be  proportional  to  the  temperature  adjustment, 

?*AA*  =  -A  Tk  (5a) 

with  a  single  proportionality, 

max(r^)  -  min(7’fc) 

T‘~1'“max(5,)-niin(5,r  <5b) 

During  the  adjustment,  the  property  of  the  water  mass 
changes  if  the  T-S  relation  does  not  preserve,  and  the 
horizontal  kinetic  energy  changes  if  the  potential  energy 
does  not  conserve.  All  of  these  changes  are  not  real. 

Thus,  it  is  urgent  to  design  a  scheme  that  satisfies  the 
three  requirements:  1)  conservation  of  heat,  salt,  and 
potential  energy;  2)  minimal  adjustment;  and  3)  pres¬ 
ervation  of  ( T ,  A)  coherency.  A  variational  adjustment 
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(VA)  scheme  with  satisfaction  of  all  the  three  require¬ 
ments  is  developed  in  this  study.  This  scheme  shows  much 
better  performance  than  the  existing  CA,  MA,  and  CMA 
schemes.  The  rest  of  the  paper  is  outlined  as  follows. 
Section  2  describes  the  VA  scheme.  Section  3  introduces 
the  data  used  in  this  paper.  Section  4  presents  a  com¬ 
parison  among  the  CA,  MA,  CMA,  and  VA  schemes. 
Section  5  shows  the  evolution  of  an  upper-ocean  mixed 
layer  with  the  four  different  schemes.  Section  6  presents 
the  conclusions. 


and  preservation  of  the  T-S  relationship.  Here,  H  is 
the  maximum  vertical  extension  of  the  profile  data; 
(A  T,  A  .S')  and 

Ap  =  /3AS  —  aAT  (8) 


are  the  adjustments  of  (7’,  S,  p)  at  z.  The  conservation 
expressions  in  Eq.  (7)  can  be  discretized  by 


V(A7*  +  A7,+1)  _ 

A,  9  V*  z, 

k=  1  z 


k+ 1 


)  =  0, 


(9a) 


2.  VA  scheme 

Without  knowing  or  calculating  the  error  covariance 
matrix  C,  minimization  of  the  simple  cost  function 

J  =  -AxTIAx  (6) 

2  v  ’ 


‘g(AS,.  +  AW(  ),Q,  and 

k=\  L 


V  (A  Pk  +  Apfc+t)  ,  x=0 

^  9  zk\zk  zk+l>  u’ 

k= 1  z 


(9b) 

(9c) 


is  used.  If  the  error  covariance  matrix  C  is  given,  the 
identity  matrix  I  is  replaced  by  C.  The  minimization  is 
conducted  with  the  strong  constrains  of  the  static  stability 
in  Eq.  (2),  conservation  of  heat,  salt,  and  potential  energy, 


■o 

AT  dz  =  0, 


-H 


■0 

A5  dz  =  0, 


-H 


■0 

gzApdz  =  0, 

-H 


(7) 


which  are  rearranged  into 


K  K 

X  ai tATfr  =  0,  X  ak^Sk  =  0, 

k= 1  k= 1 


K  K 

Ug  X  bkATk  -  Pg  X  bkASk  =  0, 

k= 1  k= 1 


with 


(10) 


(Zj  -  z2)  (z1  -  z3)  (z2  -  z4) 

al= - 9 - ,a2= - 9 - ,«3= - 9 - , 


_  (zk- 2  zk)  „  _  (zk-l  zk> 
'ak- 1  o  ’uk 


U  _z1(z1-z2)  u  _Z2(Z1-Z3)  u  _z3(z2~z4)  u  _ 
Di  2  ’°2  2  ’°3  2  ’  ' 


2  ’  k  2 

_  zk-l(zk-2  —  z*)  u  —Zk(zk-\~zk) 


-’bk  = 


(11) 


The  T-S  relation  is  obtained  using  locally  weighted 
linear  regression  from  N  surrounding  casts  (Tk\Sk\j  = 


1,2,  . . . ,  N)  (Fox  et  al.  2002): 


Sk  Sk  +  ak(Tk 


Tk), 


(12) 


N 


X  bji&k  -  Sk)(T(kj)  -  rk) 

J=  l 


N 


X  hjCrfP  ~  Tk)2 

7=1 


N 


X  bjS[ 
s  -J=: 

*k  N 


(7) 


Xft, 


7=1 


N 


X  bjtk 

T  ~M 
1  k  N 


17) 


(13) 


where  Tk  and  Sk  are  the  weighted  averages  at  the  cur¬ 
rent  location  where  the  ( 7’,  .S’)  cast  is  adjusted,  {uk\  are 


the  regression  coefficients  calculated  from  the  TV  sur¬ 
rounding  casts,  and  {bj}  are  the  weights, 


b .  =  exp 


(Xj  xq) 


C yj  ~  y  o) 

T  2 


2l 


(14) 


where  (x0,  yo)  and  (x;-,  yy)  are  the  east-west  and  north- 
south  positions  of  the  current  and  the  / th  surrounding 
casts,  respectively,  and  (Lx,  Ly)  are  the  decorrelation 
scales  in  the  (x,  y)  direction  (see  Fox  et  al.  2002). 
Equation  (12)  can  be  rewritten  as 

AA/c  _  ak^Tk  =  ^k  +  ak(Tk  ~  Tk>  ~  Sk  ’  (I5) 

which  provides  a  salinity  adjustment  with  a  preserving 
T-S  relation.  Equation  (15)  shows  the  major  differences 
between  the  existing  (CA,  MA,  CMA)  schemes  and 
the  VA  scheme.  The  existing  schemes  are  fully  local 
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TABLE  1.  The  WOA98  profiles  for  the  53.5°S,  171.5°E  grid  box  before  stabilization  (from  Locarnini  et  al.  2006,  their  Table  Bl). 


k 

Depth  (m) 

T  (°C) 

S  (ppt) 

p(Sk+1,Tk+1,  Zk)  (kg  m  3) 

p{SkJk ,  Zk)  (kg  m  3) 

Ek  (kg  m  3) 

i 

0 

7.1667 

34.4243 

26.9476 

26.9423 

0.0054 

2 

10 

7.1489 

34.4278 

26.8982 

26.9939 

-0.0957* 

3 

20 

7.0465 

34.2880 

26.9529 

26.9443 

0.0085 

4 

30 

7.0050 

34.2914 

27.0104 

26.9990 

0.0114 

5 

50 

6.9686 

34.2991 

27.0967 

27.1028 

-0.0061* 

6 

75 

7.0604 

34.3073 

27.2406 

27.2120 

0.0286 

7 

100 

6.9753 

34.3280 

27.3892 

27.3560 

0.0332 

8 

125 

6.9218 

34.3604 

27.5164 

27.5046 

0.0117 

9 

150 

6.8919 

34.3697 

27.6000 

27.6316 

-0.0316* 

10 

200 

6.9363 

34.3364 

27.8123 

27.8302 

-0.0179* 

11 

250 

7.0962 

34.3415 

28.0295 

28.0421 

-0.0126* 

12 

300 

7.1622 

34.3367 

28.2684 

28.2593 

0.0092 

13 

400 

6.8275 

34.2852 

28.6664 

28.7281 

-0.0618* 

14 

500 

7.4001 

34.3123 

29.3699 

29.1238 

0.2461 

15 

600 

6.2133 

34.4022 

29.9386 

29.8292 

0.1094 

16 

700 

5.9186 

34.4868 

30.5869 

30.3978 

0.1891 

17 

800 

4.5426 

34.4904 

31.0754 

31.0488 

0.0266 

18 

900 

4.1263 

34.4558 

31.6539 

31.5377 

0.1162 

19 

1000 

3.3112 

34.4755 

32.1176 

*  Static  instability. 


adjustments;  that  is,  no  surrounding  ( T ,  S )  casts  are  used. 
However,  the  VA  scheme  uses  N  surrounding  ( T ,  S ) 
casts  to  preserve  the  T-S  coherency. 

A  successive  quadratic  programming  method 
(Schittkowski  1986)  is  used  to  solve  the  general  nonlinear 
optimization  problem.  When  optimality  is  not  achieved,  a 
positive  definite  approximation  of  the  Hessian  is  updated 
according  to  the  modified  Broyden-Fletcher-Goldfarb- 
Shanno  method  (Powell  1978). 

3.  Data 


temperature  data  than  salinity  data  on  the  New  Zealand 
Plateau  for  October.  This  contributes  to  six  small  (on  the 
order  of  10  2  kg  m  3 )  inversions  in  the  local  potential 
density  field  calculated  from  objectively  analyzed  tem¬ 
perature  and  salinity  fields  (Table  1).  This  ( T ,  S )  cast  is 
used  for  evaluating  the  performance  of  the  VA  scheme 
versus  the  existing  CA,  MA,  and  CMA  schemes. 

4.  Verification 

a.  Metrics  for  verification 


Evident  false  static  instability  was  found  in  earlier 
versions  of  the  World  Ocean  Atlas  (WOA)  such  as  more 
than  44%  in  1982  version  of  the  WOA  ( WOA82 ;  Jackett 
and  McDougall  1995)  and  around  3%  in  WOA94  (Chu 
and  Fan  2010a).  The  static  instability  disappears  in  re¬ 
cent  versions  of  WOA  due  to  the  utilization  of  the  mod¬ 
ified  MA  scheme,  which  adjusts  the  ( T ,  .S’)  profiles.  An 
example  as  described  in  appendix  B  of  Locarnini  et  al. 
(2006)  is  used  for  illustrating  the  capability  of  the  scheme. 
The  area  chosen  for  this  example  is  the  1°  latitude- 
longitude  box  centered  at  53.5°S,  171. 5°E  from  the  1998 
version  of  the  World  Ocean  Atlas  ( WOA98 ).  The  area  of 
interest  is  on  the  New  Zealand  Plateau,  with  a  bottom 
depth  below  1000  m  and  above  1100  m.  The  month  is 
October,  during  the  early  austral  summer.  There  is  no 
temperature  or  salinity  data  within  the  chosen  1°  box. 
Thus,  the  objectively  analyzed  values  in  this  1°  box  will 
be  dependent  on  the  seasonal  objectively  analyzed  field 
and  the  data  in  nearby  1°  grid  boxes.  There  is  much  more 


Several  metrics  are  used  to  evaluate  the  magnitude  of 
the  (T,  S )  adjustment  and  the  conservation  of  heat,  salt, 
and  potential  energy.  Similar  to  Chu  and  Fan  (2010b), 
the  relative  root-mean  adjustment  (RRMA)  is  used  to 
identify  the  changes  in  temperature, 


RRMAj, 


I  (A  Tkf 
k= 1 


max(rfc)-min(rfc)’ 


(16a) 


salinity, 


RRV1A, 


I  (AS,)2 

k=l 


max(S,)  -  min(S,)’ 


(16b) 


and  temperature-salinity, 
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Temp  difference:  0.0000  Temp  difference:  0.0000 


Fig.  1.  Original  (dashed)  and  adjusted  (solid)  temperature  profiles  Tk  at  the  53.5°S,  171.5°E  grid 
box  using  the  (a)  VA,  (b)  CA,  (c)  MA,  and  (d)  CM  A  schemes. 
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Table  2.  RRMAj  and  RRMAj  with  the  four  adjust  schemes  on 
the  original  profile  shown  in  Table  1.  Here,  the  values  for  CA,  MA, 
and  CMA  are  from  Chu  and  Fan  (2010a). 

Schemes 

RRMAr 

RRMAj 

CA 

0.0287 

0.1905 

MA 

0.0220 

0.0492 

CMA 

0.0241 

0.0241 

VA 

0.0043 

0.1192 

RRMA  = 


ma x(Tk)  -  min (Tk)  max(Sk)  —  min(5A.) ' 

(16c) 


is  used  as  a  metrics  to  identify  the  change  of  the  adjusted 
from  original  casts.  RRMA  represents  the  mean  ad¬ 
justment  relative  to  the  range  of  a  profile.  The  total  heat, 
salt,  and  potential  energy  (PE)  changes  of  the  water 
column  within  this  1°  X  1°  grid  box  are  estimated  by 


A  Q=APocp 


•o 

A  Tdz, 


-H 


APE  = 


■o 

gzAp  dz , 
-H 


A(Salt)  =  A 


■o 


A  Sdz, 

-H 


(17) 


where  p0  (=1028  kg  m~3)  is  the  characteristic  density, 
cp  (=4002  J  kg-1  KT1)  is  the  specific  heat  for  the  sea¬ 
water,  H  =  1000  m,  and  A  is  the  area  of  the  grid  box, 


Here,  R  (=6370  km)  is  the  earth  radius,  and  cp  (=53.5°) 
is  the  latitude  of  the  grid  box. 

b.  Verification 

Figure  1  shows  the  comparison  of  temperature  ad¬ 
justments  using  the  four  schemes  (CA,  MA,  CMA,  and 
VA).  The  VA,  CA,  and  CMA  schemes  conserve  the 
heat,  but  the  MA  scheme  does  not.  The  three  schemes 
with  a  minimum  adjustment  requirement  (VA,  MA,  and 
CMA)  generate  the  new  profiles  by  basically  following 
the  original  profile  pattern,  but  the  CA  scheme  does 
not.  Among  the  three  heat  conserved  schemes,  the  CA 
scheme  has  the  largest  change  (0.0296)  and  the  VA 
scheme  has  the  smallest  change  (0.0043)  (see  Table  2). 
For  the  layer  between  400  and  600  m,  the  MA  scheme 
has  largest  adjustment  among  VA,  CA,  MA,  and  CMA, 
which  may  lead  to  the  largest  RRMAr  while  the  other 
two  schemes  (VA  and  CMA)  are  nearly  consistent  with 
the  original  profile. 


Figure  2  shows  the  comparison  of  salinity  adjustment 
using  the  four  schemes  (CA,  MA,  CMA,  and  VA).  The 
VA,  CA,  and  CMA  schemes  conserve  the  salt,  but  the 
MA  scheme  does  not.  The  three  schemes  with  a  mini¬ 
mum  adjustment  requirement  (VA,  MA,  and  CMA) 
generate  the  new  profiles  by  basically  following  the  origi¬ 
nal  profile  pattern,  but  the  CA  scheme  does  not.  Among 
the  three  salt  conserved  schemes,  the  CA  scheme  has  the 
largest  change  (0.1905),  and  the  CMA  scheme  has  the 
smallest  change  (0.0241),  while  the  VA  scheme  has 
a  modest  change  (0.1192)  (see  Table  2).  Figure  3  shows 
the  comparison  of  the  static  stability  (i.e.,  N2)  adjust¬ 
ment  using  the  four  schemes  (CA,  MA,  CMA,  and  VA). 
Among  them,  the  CMA  and  MA  schemes  have  small 
adjustments,  and  the  VA  and  CA  schemes  have  rela¬ 
tively  large  adjustments.  It  is  also  noted  that  salinity  for 
the  VA  method  is  adjusted  over  the  entire  profile  and 
the  static  stability  increases  at  800-m  depth  where  the 
buoyancy  frequency  is  near  neutral  (Fig.  3a).  However, 
the  salinity  for  the  existing  (CA,  MA,  and  CMA)  methods 
is  adjusted  above  600-m  depth.  This  may  be  caused  by 
the  use  of  N  surrounding  casts  in  the  VA  scheme  [see 
Eq.  (15)]. 

Table  3  lists  the  overall  features  of  the  adjustment 
using  the  four  different  schemes.  The  RRMA  using  VA 
(0.1235)  is  smaller  than  that  using  CA  (0.2192),  and 
larger  than  that  using  either  MA  (0.0712)  or  CMA 
(0.0482).  The  heat,  salt,  and  potential  energy  are  fully 
conserved  using  the  VA  scheme,  but  are  not  fully  con¬ 
served  using  the  other  schemes.  It  causes  a  heat  loss  of 
1.4  X  108  J  m~2,  a  salt  loss  of  2.31  kg  m~z,  and  a  po¬ 
tential  energy  gain  of  30  700  J  m  2  using  the  MA  scheme, 
plus  a  potential  energy  gain  of  4600  J  m  2  using  the 
CMA  scheme,  and  a  potential  energy  gain  of  32  800  J  m  2 
using  the  CA  scheme.  With  the  geostrophic  adjustment, 
the  potential  energy  can  been  consumed  by  either  the 
production  of  kinetic  energy  or  gravitational  waves  that 
travel  to  infinity,  effectively  carrying  energy  away  from 
the  region  (Marotzke  and  Scott  1999).  Thus,  it  is  likely 
that  the  nonconservation  of  the  potential  energy  during 
the  adjustment  may  change  the  velocity  field. 


5.  Free  convection 

a.  Theoretical  limit 

Let  us  now  consider  free  convection  (no  winds)  in 
a  resting  ocean  of  constant  stratification  V,h  (subscript 
“th”  for  thermocline)  with  a  uniform  and  widespread 
buoyancy  loss  ( B0 )  at  the  surface.  The  response  to  wide¬ 
spread  cooling  is  one  in  which  relatively  small  convection 
cells  (plumes)  develop.  Fluid  parcels  in  contact  with  the 
surface  will  become  dense  and  sink  under  gravity,  driving 
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Fig.  2.  Original  (dashed)  and  adjusted  (solid)  salinity  profiles  Sk  at  the  53.5°S.  171. 5°E  grid  box 
using  the  (a)  VA,  (b)  CA,  (c)  MA,  and  (d)  CMA  schemes. 
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Fig.  3.  Original  (dashed)  and  adjusted  (solid)  N 2  profiles  at  the  53.5°S,  171.5°E  grid  box  using  the 
(a)  VA,  (b)  CA,  (c)  MA,  and  (d)  CMA  schemes. 
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Table  3.  RRMA  and  changes  in  heat,  salt,  and  potential  energy 
using  the  four  adjust  schemes  on  the  original  profile  shown  in 
Table  1. 


Scheme 

RRMA 

Heat  change 
per  unit  area 
(10s  J  nV2) 

Salt  change 
per  unit  area 
(kg  m"2) 

Potential  energy 
change  per  unit 
area  (J  nV2) 

CA 

0.2192 

0 

0 

32  800 

MA 

0.0712 

-1.4 

-2.31 

30  700 

CMA 

0.0482 

0.0 

0.0 

4600 

VA 

0.1235 

0.0 

0.0 

0.0 

the  “free  convective  layer”  below.  Buoyancy  is  drawn 
upward,  across  the  convective  layer,  offsetting  its  loss 
from  the  surface  (Marshall  and  Schott  1999).  For  the 
same  buoyancy  loss  at  the  ocean  surface,  the  free  con¬ 
vection  depends  on  the  temperature  and  salinity  pro¬ 
files,  which  are  adjusted  with  various  schemes  (i.e.,  CA, 
MA,  CMA,  and  VA).  Thus,  it  is  reasonable  to  use  this 
process  to  investigate  the  impact  of  various  adjustment 
schemes  on  the  evolution  of  the  free  convective  layer. 
If  entrainment  of  stratified  fluid  from  the  base  of  the 
mixed  layer  can  be  neglected,  the  time  evolution  of  the 
mixed  layer  depth  (MLD)  is  calculated  by  (Turner 
1973) 


M0  =  tt\  2 


N, 


th 


Bq  dt. 


(18) 


If  Bq  is  constant  in  time,  Eq.  (18)  becomes 


mo=t rv2Bo{- 

Jvth 


(19) 


In  real  oceans,  the  stratification  below  the  mixed  layer  is 
not  constant  and  the  entrainment  of  stratified  fluid  from 
the  base  of  the  mixed  layer  cannot  be  neglected.  The 
time  evolution  of  MLD  does  not  follow  (18)  or  (19). 
Thus,  numerical  modeling  is  needed. 

b.  Mixed  layer  modeling 

A  one-dimensional  mixed  layer  model  (Price  et  al. 
1994)  is  used  to  predict  the  evolution  of  the  ocean  mixed 
layer  with  free  convection  in  a  horizontally  homogenous 
ocean.  In  this  model,  turbulent  mixing  occurrence  in¬ 
volves  three  stability  criteria.  The  vertical  mixing  is  ac¬ 
counted  for  by  static  instability,  and  bulk  and  gradient 
Richardson  number  mixing.  The  critical  bulk  and  gra¬ 
dient  Richardson  numbers  are  set  to  0.65  and  0.25,  re¬ 
spectively.  The  model  output  includes  the  time  evolution 
of  the  temperature  and  salinity  profiles  [T(z,  t),  S(z,  /)]  as 
well  as  MLD  [fi(f)]. 


c.  Impact  of  the  schemes 

The  five  ('/’,  S)  profile  pairs  (one  original  and  four 
adjusted,  shown  in  Figs.  1  and  2)  do  not  have  constant 
stratification  in  the  thermocline.  So,  the  theoretical  limit 
in  (18)  or  (19)  might  not  be  suitable  for  predicting  h(t). 
The  mixed  layer  model  (Price  et  al.  1994)  is  used  to 
simulate  free  convection  due  to  a  steady  heat  loss  of 
220  W  m  2  at  the  sea  surface.  The  model  is  integrated 
for  5  days  with  the  original  and  adjusted  temperature 
and  salinity  profiles  (by  the  four  schemes)  as  the  initial 
conditions. 

Significant  differences  are  found  in  the  model- 
simulated  5-day  time  evolution  patterns  of  vertical  tem¬ 
perature  (Fig.  4)  for  the  original  and  the  four  adjusted 
(CA,  MA,  CMA,  and  VA)  profiles  as  the  initial  condi¬ 
tions.  Hereafter,  the  adjusted  profiles  will  be  called  the 
CA,  MA,  CMA,  or  VA  profiles  for  convenience.  The 
temperature  field  is  similar  in  the  mixed  layer  between 
using  the  original  and  MA  profiles  up  to  day  3.5  and 
different  afterward  with  a  vertically  uniform  layer  (down 
to  500-m  depth)  for  the  MA  profile  and  still  stratified 
below  100  m  for  the  original  profile.  The  temperature 
field  below  the  mixed  layer  is  nearly  0.5°C  colder  for  the 
MA  profile  than  that  for  the  original  profile  within  about 
3.5  days.  At  day  3.5,  the  thermocline  suddenly  breaks 
down  for  the  MA  profile,  which  contributes  to  the  great 
mixed  layer  deepening.  The  temperature  is  similar  in  the 
mixed  layer  and  there  is  little  difference  (0.1°-0.2°C) 
below  the  mixed  layer  between  using  the  original  and 
CMA  profiles  for  the  whole  5  days  of  simulation.  The 
temperature  is  slightly  warmer  in  the  mixed  layer  using 
the  VA  profile  than  the  original  profile,  but  it  is  very 
similar  below  the  mixed  layer  for  the  two  profiles.  Com¬ 
pared  with  the  other  three  schemes  and  the  original  pro¬ 
file,  using  the  CA  profile  reveals  significant  differences. 
The  temperature  in  the  upper  500  m  throughout  the  re¬ 
gion  remains  homogeneous  without  marked  stratification. 

The  simulated  time  evolution  patterns  of  MLD,  h(t), 
are  very  different  among  the  original  and  four  adjusted 
initial  profiles  (Fig.  5).  The  MLD  increases  linearly  with 
almost  the  same  rate  between  the  CMA  and  the  original 
profiles,  although  the  MLD  is  25  m  thicker  for  the  CMA 
profile.  The  MLD  increases  with  time  similarly  when 
comparing  the  MA  and  VA  profiles  within  day  1  and 
then  very  differently  afterward.  The  MLD  suddenly  in¬ 
creases  from  100  to  500  m  at  day  3.5  for  the  MA  profile; 
however,  the  MLD  increases  to  almost  the  same  level  as 
the  theoretical  limit  (19)  with  a  given  Vt h  (  =  1.414  X 
10~3  s_1)  for  the  VA  profile.  This  may  indicate  a  posi¬ 
tive  feature  of  the  VA  scheme.  For  the  CA  profile,  the 
initial  MLD  is  deeper  with  about  500  m,  which  may  in¬ 
hibit  the  mixing  development  due  to  the  cooling.  During 
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Fig.  5.  The  variability  of  mixed  layer  depth  during  5  days  of 
simulation  for  the  original  profile  conditions  and  four  adjustment 
schemes  of  static  instability,  as  well  as  theory  prediction. 


Fig.  6.  The  variability  of  gravitational  potential  energy  per  unit 
area  relative  to  initial  state  during  5  days  of  simulation  for  the 
original  profile  conditions  and  four  adjustment  schemes  of  static 
instability  (J  m~2). 


5-day  simulations,  the  MLD  represents  the  evident  os¬ 
cillation  of  the  1-day  period. 

As  described  by  Huang  and  Wang  (2003),  with  a  lim¬ 
ited  amount  of  heat  loss  due  to  cooling,  a  thin  layer  of 
water  at  the  top  of  the  mixed  layer  becomes  heavier  than 
the  water  below.  This  unstable  stratification  leads  to  a 
convective  adjustment  when  gravitational  potential  en¬ 
ergy  is  transformed  into  the  small-scale  energy,  most  of 
which  is  eventually  dissipated  into  thermal  energy.  This 
energy  loss  is  probably  one  of  the  most  important  compo¬ 
nents  of  the  energetic  balance  in  world  oceans.  Figure  6 
shows  the  simulated  time  evolution  of  the  potential  energy 
per  unit  area  relative  to  its  initial  value  for  the  original  and 
four  adjusted  profiles.  For  the  original  profile,  it  suddenly 
decreases  about  500  J  m  2  at  the  very  beginning  of  the 


cycle  due  to  the  release  of  the  static  instability  at  the  initial 
stage,  and  then  decreases  linearly  to  1800  J  m  2  at  the  end 
of  day  5.  For  the  MA  profile,  it  decreases  gradually  within 
about  3.5  days,  but  very  rapidly  afterward  increases  to 
3300  J  mT2  on  day  4.  For  both  the  CMA  and  VA  profiles,  it 
decreases  gradually  due  to  the  surface  cooling  during  the 
5  days  of  simulation.  After  5  days  of  free  convection  sim¬ 
ulation,  the  potential  energy  reduction  is  1615  J  in  2  for  the 
VA  profile  (minimum),  1739  J  in  2  for  the  CMA  profile, 
1883  J  mT2  for  the  original  profile,  4500  J  in  2  for  the  MA 
profile,  and  6410  J  tn  2  for  CA  profile  (maximum). 

6.  Conclusions 

A  new  fully  conserved  minimal  adjustment  scheme  is 
developed  to  eliminate  the  static  instability  of  raw  and 
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averaged  observational  hydrographic  data.  This  method 
adjusts  the  temperature  and  salinity  profiles  {A 7^,  AS*,  k  = 
1,2,  .. .  ,K}  simultaneously  and  efficiently  on  the  base  of 
three  types  of  constraints:  (a)  heat,  salt,  and  potential 
energy  conservation;  (b)  a  preserved  T-S  relationship; 
and  (c)  the  removal  of  static  instability  by  minimal  ad¬ 
justment  of  the  T-S  casts.  With  these  constraints,  a  var¬ 
iational  algorithm  is  designed  to  remove  the  false  static 
instability  induced  by  temperature-salinity  data  analysis/ 
assimilation. 

The  VA  scheme  can  contribute  to  the  data  quality 
control  process  in  data  assimilations  since  it  does  not 
simply  reject  profiles  with  static  instability.  This  method 
edits  the  profiles  with  the  inequality  constraint  to  satisfy 
static  stability.  Such  a  quality  control  process  yields 
a  reasonable  response  to  a  one-dimensional  mixing 
model. 

In  principle,  variational  data  assimilation  can  require 
the  constraint  conditions  to  satisfy  static  stability.  These 
constraints  suggested  in  the  present  paper  can  be  forced 
directly  to  the  cost  function  of  the  variational  data  as¬ 
similation  in  order  to  gain  the  assimilation  result  of  static 
stability.  The  quadratic  norm  AxTCAx  can  be  used  to 
construct  a  rejection-acceptance  criterion.  With  the 
features  of  (a)  full  (heat,  salt,  and  potential  energy) 
conservation,  (b)  minimal  adjustment,  and  (c)  ( T ,  S ) 
coherency,  the  VA  scheme  can  be  used  in  global  ocean 
(T,  S )  data  analysis  and  assimilation.  However,  further 
examination  is  needed  since  the  results  may  vary  with 
changing  norms  and  since  only  limited  data  are  used  in 
this  study  for  the  verification. 
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